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nteractive Physics (IP) provides the 

user with a computer-based labora- 

tory that simulates the mechanical 
behavior of objects in a two-dimension- 
al world. Its most seductive feature is 
that the primary method of data input 
and output is graphical. For the most 
part, the user literally draws the experi- 
ment and its initial conditions. JP then 
calculates and displays the consequent 
action in the form of an animation. 

Providing such capability is a tall 
order indeed—at least if one cares to get 
the physics right—and JP is not entirely 
successful. Nevertheless, its usefulness 
as a demonstration aid and as an 
instrument for a myriad of student 
assignments and mini-research projects 
make it well worth the cost. JP will 
teach anyone some physics. 

In this review, I intend to present a 
reasonably detailed examination of how 
IP works in order to expose both its 
strengths and weaknesses. To prepare 
readers who are unfamiliar with JP for 
what follows, the remainder of this 
section is devoted to an overview of the 
program’s operation. Readers with pre- 
vious experience of using JP may prefer 
to skip to the next section. 

In preparation for drawing an 
experiment, the user adjusts properties 
of the “world” in which it is to take 
place. The world’s size may be adjusted 
from as large as one square kilometer 
down to one-hundredth of a square 
meter and scroll bars and zoom com- 
mands can be used to view all or any 
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part of the world. A “downward” (i.e. 
toward the bottom of the screen) 
gravitational field vector can be adjust- 
ed in magnitude from 0 to 25 N/kg, and 
a very rudimentary form of quadratic 
air resistance can be invoked by setting 
what amounts to a universal drag 
coefficient for which no units are 
specified. (See the discussion of air 
resistance in the section entitled “Could 
It Be Any Better?’’.) JP provides 
toggles for the display of coordinate 
axes, a coordinate grid, and horizontal 
and vertical rulers. 

The user creates an experiment by 
drawing “mass objects” in the form of 
squares, rectangles, circles and poly- 
gons, using standard Macintosh draw- 
ing tools. A “grid snap” option is 
provided, allowing for more control 
over the positions and sizes of objects by 
forcing their vertices (horizontal and 
vertical extrema in the case of circles) 
to occupy points on a coarse lattice. 
Once they’re drawn, the objects can be 
resized, repositioned, rotated and co- 
lored. Each object carries a user-specifi- 
able name, and is characterized by an 
adjustable mass and two adjustable— 
but unquantified—parameters, which 
IP calls the “Surface Parameter” and 
the “Elasticity.” 

Whenever two objects come into 
contact, these parameters work in com- 
bination to determine, in an empirical 
manner, the frictional and rebound 
characteristics of the interaction. Alter- 
natively, one may select specific pairs of 
objects and explicitly specify their mu- 
tual coefficients of friction. The initial 
center-of-mass position, rotational ori- 
entation, and translational and rota- 
tional velocities of each mass object can 
be set in a spec sheet accessed by 
double-clicking on the object. When 
precision is less important, the initial 
translational velocity of an object may 
be specified graphically by selecting the 
object and dragging (from its center-of- 
mass) a vector of the desired size and 
direction. Any object may be anchored 
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in space to form a fixed platform or 
barrier. ; 

IP can display vectors representing 
an object’s velocity or acceleration; the 
individual forces exerted on it by other 
objects, air resistance or gravity; or the 
net force acting on it. In addition, the 
user can request meters which monitor 
these quantities as well as time, posi- 
tion, rotational position, velocity or 
acceleration, translational or rotational 
kinetic energy, translational or rota- 
tional momentum, net torque, and the 
forces of contact between any two 
selected objects. The meters can pro- 
vide digital or analog readouts or strip- 
chart-type recordings of the monitored 
value. 

The user may use drawing tools to 
connect mass objects to each other 
and/or to fixed points in space with 
“no-mass objects”—strings, springs, 
and dashpots. These objects have ap- 
propriate adjustable parameters, which, 
again, can be accessed by double- 
clicking or via a menu item. They apply 
their forces at the points to which their 
ends are connected. No-mass objects do 
not interact with objects other than 
those to which they are connected, so 
that, for instance, one cannot form a 
pulley by wrapping a string around a 
circle. 

A final drawing tool produces 
arbitrary constant force vectors which 
can be applied at any point on mass 
objects and which can maintain either 
their absolute orientation in space, or 
their orientation with respect to the 
object itself. 

When finished with the set-up of 
an experiment, the user clicks on 
“Run”, and JP animates all objects in 
accordance with Newton’s laws, using a 
simple fixed time-step numerical inte- 
gration method. (See the next section 
for a brief discussion of the computa- 
tional algorithm.) The user can specify 
the integration time-step. As with all 
numerical integration schemes, the sim- 
ulation accuracy (generally) increases 
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with smaller time-steps at the necessary 
expense of simulation speed. 

Because substantial amounts of 
time are spent doing the physics, ZP 
records animation frames in a “tape 
player memory” which can be used to 
replay the animation forward or back- 
ward, either continuously (at increased 
speed) or one step at a time. Tape 
player memory can be allocated in 
chunks from 25 kbytes up to the limit 
imposed by the amount of memory 
available to the program. The default 
allocation of 400 kbytes is sufficient to 
store about 4000 frames of a simple, 
two-object experiment, but will manage 
only about 120 frames of the most 
complex experiments. An optional 
warning dialog box pops up when the 
memory is first filled and, if the user 
allows JP to proceed, new frames will 
overwrite the earliest frames, preserv- 
ing only the most recent segment of the 
experiment. A “Skip Frames” option 
allows the user to view only one frame 
out of every two, four, eight or 16 to 
further speed up replayed experiments. 
Finally, a “Tracking” option leaves 
unerased every first, second, fourth, 
eighth, sixteenth or thirty-second 
frame, giving a stroboscopic picture of 
the experiment. Any frame (or tracked 
sequence of frames) can be printed. 

At the end of an experiment, the 
user may “export” the data collected by 
all meters, to a data document which 
can be analyzed later using a spread- 
sheet or some other type of data 
analysis program. Alternatively, data 
from any set of the available meters 
may be copied to the clipboard and 
immediately “pasted” into the analysis 
program. 

IP comes with two folders of 
example scenarios. The “Demonstra- 
tion” folder contains a set of simple 
experiments providing a quick tour of 
IP’s capabilities. The “Experiments” 
folder offers experiments in one- and 
two-dimensional kinematics, dynamics, 
energy, momentum, collisions, rota- 
tional dynamics, equilibrium, and oscil- 
lations which are intended as the basis 
for a variety of student “lab” exercises. 
These scenarios stick to relatively sim- 
ple situations and are probably most 
useful as jumping-off points for student 
investigations or demonstrations. A 
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supplementary book entitled Physics 
Interactions serves as a lab manual for a 
set of experiments on an accompanying 
disk. 


How Does It Do It? 


As is true for any numerical simulation, 
the user must be alert for unphysical 
artifacts in the output of JP that result 
from the fundamental limitations of its 
computational algorithms. Without 
some understanding of what those 
algorithms are, one can’t make fully 
intelligent use of JP and will inevitably 
be frustrated by a number of ill- 
mannered and unphysical behaviors. 

For its routine operations, JP uses 
what amounts to a modified, fixed time- 
step, Eulerian integration scheme. At 
the beginning of each integration step, 
IP determines the individual forces and 
their points of application, and uses this 
information to find the net force and net 
torque acting on each object. Dividing 
these quantities by the mass and rota- 
tional inertia gives the translational and 
rotational accelerations. The transla- 
tional and rotational impulses are cal- 
culated—under the assumption that the 
forces and torques are constant over the 
duration of the next integration time- 
step—and added to the old transla- 
tional and rotational momenta to ob- 
tain new values. Dividing these mo- 
menta by the mass and rotational 
inertia of the object provides the new 
translational and rotational velocities. 

Finally, ZP calculates the new 
translational and rotational positions, 
assuming that the new values for the 
velocities were in effect throughout the 
previous time-step. This modification of 
the standard Euler method has been 
shown to yield more accurate results for 
a given step size in many cases of 
interest. ' 

The foregoing procedure would be 
simply executed if it were not for 
complications that arise primarily in 
determining the appropriate magni- 
tudes and effective points of application 
for forces of contact. These problems 
only begin with trying to determine 
when and where pairs of objects come 
into contact in the first place. They 
become especially acute when one cares 
to model the effects of friction and 


partial elasticity with any semblance of 
physical correctness. However, in many 
circumstances, the last straw is the 
commitment to do all of this with a 
constant integration step size for the 
sake of both maximizing and maintain- 
ing animation speed. Forces of contact 
are particularly unforgiving in this 
regard. 

IP adopts several empirical rules 
which, in effect, assert desired out- 
comes to determine requisite forces. 
For instance, at least in the simple case 
of a collision between two objects, the 
product of the elasticity parameters 
determines the ratio of the relative 
velocity after the collision to that before 
the collision. During the first integra- 
tion step after ZP detects a collision— 
which will, in general, be after the 
drawn boundaries of the two objects 
visibly overlap—it applies equal and 
opposite normal forces calculated to 
provide the impulse necessary to yield 
the desired final velocities in a single 
time-step. 

Additional complications arise 
from JP’s need to remember that it has 
“taken care of” the collision even 
though the boundaries will, in general, 
continue to overlap for one or more 
frames. Of course, it can simply deter- 
mine whether the relative motion of the 
objects is toward or away from each 
other, but what happens if other objects 
are simultaneously interacting—push- 
ing each other together? As you may 
see, it can get complicated. 

Another empirical rule is used to 
treat initially-stretched ropes. When a 
string is found to be longer than its 
resting length, ZP applies a single-step 
impulse to whatever objects may be 
found connected to its ends. The im- 
pulse is calculated on the basis of a 
desire to return the string to its resting 
length in a relatively small number of 
frames and, therefore, depends not only 
on how much the string is stretched, but 
also on the masses and velocities of the 
objects to which it is attached, and the 
integration time step. Furthermore, if 
IP determines that the string is making 
sufficiently rapid progress toward its 
resting length, the string will apply no 
further forces to the objects, regardless 
of the fact that it may still be substan- 
tially stretched. In situations where 


other forces work against the desires of 
the string, ZP will reach a compromise 
in which the string exerts whatever 
force is necessary to effect a damped 
return to its resting length. Although 
one may find this algorithm disturbing 
on physical grounds, in most cases it is 
probably preferable to a more realistic 
one in which the string would merely 
break if it were stretched by any 
noticeable amount. 


Any Problems? 


In the course of putting ZP through its 
paces, I discovered a number of bugs. In 
a conversation with David Baszucki, 
president of Knowledge Revolution 
and author of Interactive Physics, 1 was 
informed that the most important four 
of these have been exterminated and 
will not be found in Version 1.2.1, 
which is now shipping. He also told me 
that registered owners of Version 1.2 
will receive the update free of charge. 
The four exterminated bugs were char- 
acterized by the following pernicious 
behaviors: 

e Causing the effective coefficient of 
friction between two objects to behave 
as though it were four times what the 
user had specified. 

è Losing track of the area of objects 
which have been resized. This error also 
affected the calculation of the object’s 
rotational inertia. 

e Including two spaces before each 
numerical datum (only one space in 
front of negative values) in the tab- 
delimited data export files. Although 
not really a bug, these spaces caused 
some spreadsheets and data analysis 
programs to consider the data as text, 
and necessitated some form of conver- 
sion before numerical manipulations 
could be accomplished. 

e Failures of the contact detection 
algorithm that occasionally allowed 
objects under certain seemingly benign 
circumstances to pass directly through 
each other. 

Other bugs I encountered—many 
of which the people at Knowledge 
Revolution were apparently previously 
unaware of—included the following: 
ə In at least one rather special situation, 
the frictional interaction between two 
objects acts in a direction which would 


increase the relative velocity between 
the points of contact. 

e Another failure of the contact detec- 
tion algorithm that sometimes occurs 
when a collision between two objects 
takes place at vertices of both objects. 
eè Displaying reversed normal and fric- 
tion force vectors in certain circum- 
stances. This appears to be the result of 
the program becoming confused about 
which of a given contacting pair has 
been selected for vector display. Saving 
such an experiment and reopening it 
causes the vectors to get assigned 
properly. 

e A tendency for the air force vector to 
lock up when an experiment is run 
backwards, and to disappear if an 
experiment is reset and run from the 
Start. 

eA lack of communication between the 
surface parameter that is displayed for 
individual objects, and the coefficient of 
friction that is displayed when two 
objects are simultaneously selected. 
When one sets the surface parameters 
and subsequently examines the coeffi- 
cient of friction between the two ob- 
jects, it will be displayed as zero even 
though running the experiment reveals 
that it certainly is not. 

e Reminiscing strip charts. On the 
initial run of an experiment with a strip- 
chart meter, the strip-chart will often 
lose track of the current value being 
monitored and begin replaying data 
from the beginning of the experiment. 
The problem usually disappears if the 
experiment is reset and run again. 

@ WYSIWYG failure. Printed versions 
of experiments are not always identical 
to screen versions. In one case, a vector 
labeled (correctly) “FN” on the screen 
was printed as “FG.” 

Note that most of these are rela- 
tively benign bugs that do not affect the 
simulation itself. Anyone working with 
IP for a reasonable amount of time will 
stumble onto other odd behaviors. 
Most of the time these are not so much 
bugs as they are simple artifacts of the 
numerical methods used by JP to model 
interactions. Many of these artifacts 
can be tamed simply by reducing the 
integration step-size. For instance, if 
the single-step impulse on an object due 
to air resistance is larger than its 
current momentum, then air resistance 


acting alone will instantly reverse the 
object’s direction of motion. Reducing 
the integration step size will reduce the 
very large error caused by JP’s assump- 
tion that the drag force is constant over 
that interval. 

A more subtle artifact sometimes 
emerges if the integration step is too 
small. Evidently in the interest of 
calculational speed, JP truncates most 
kinematic data at a predetermined 
decimal place intended to limit the 
number of significant digits that it has 
to carry around. As a result, if an 
impulse is non-zero, but too small to 
cause a one-digit change in the related 
momentum in a single step, then the 
momentum will not be altered by the 
impulse. Short of changing the forces 
on, or inertia of, the object, only 
increasing the integration time step will 
eliminate this pathological behavior. 
Truncation error also pops up in situa- 
tions where small consecutive changes 
in a parameter may not affect the 
experiment at all until the parameter is 
changed by a critical amount. At that 
point, a substantial alteration of the 
experiment may be observed. 

The moral of these stories is that 
one must be aware of the algorithms 
used and be alert for their effects—a 
fact of life in all computational work. 


Any Tips? 


Experiments involving objects which 
slide along each other’s surfaces are 
complicated by initial bounces that 
result from not being able to position or 
rotate the objects so that their surfaces 
are precisely aligned and just touching. 
There is, however, a nice way to obtain 
such initial conditions: turn the elasti- 
city parameters of both objects way 
down and let the experiment run just 
long enough for the bouncing to cease. 
Stop the experiment and choose “Start 
Here” from the “Experiment” menu. 
The current (bounceless) conditions 
become the new initial conditions. Now 
you can re-adjust the elasticities and 
other parameters. 

The color option may be appealing 
to users having color monitors, but 
there are few situations in which it is 
genuinely helpful, and it should be 
understood that the computational 
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overhead involved in keeping track of 
the extra color bits exacts its price in 
animation speed. In regular use, the 
monitor cdev (in the Macintosh Con- 
trol Panel) should be used to select the 


Color Mode 1 bit 


Simple task 
Complex task 





tion of this kind is essential for the 
intelligent use of JP. 

e Pivot points to provide torqueless 
rotation of an object about any point. 
(This is a necessary first step toward 


2 bit 4 bit 8 bit 


Table 1. Animation speeds (in frames per minute) for two ZP tasks running in 1-, 2-, 4- and 8-bit 
monitor modes on a Mac IIci. The “simple task” was a ball bouncing on an anchored platform. The 
“complex task” involved nine objects inside an anchored container object connected with three 


springs and a dashpot. 


1-bit (or “Black & White”) monitor 
mode. When color is required, limit 
yourself to the 2-bit (or “four-color” ) 
mode. Table 1 provides a comparison of 
approximate animation speeds on .my 
Mac IIci (in frames per minute) for 
two tasks—one simple and the other 
complex—running in the 1-bit, 2-bit, 4- 
bit and 8-bit monitor modes. Note that 
there is a much larger relative penalty 
for simple tasks than for complex 
tasks—a result of the more nearly-fixed 
time (regardless of animation complex- 
ity) required to update the display after 
each integration step. 


Could It Be Any Better? 


No question. From the moment JP was 
released, Knowledge Revolution has 
been on the receiving end of a host of 
“Why-doesn’t-it’s”, ‘‘It-needs-to’s’’, 
and “I-have-to-have’s.” A sampler of 
items from my own list follows. 
First, a number of modest requests 
(in rough order of the importance I 
would attach to them): 
e Oscillating and/or revolving drive 
points to which one can attach strings, 
springs and dashpots, and the ability to 
“track” an experiment at constant drive 
phase. (This is absolutely the number 
one request in my book. This simple 
addition would enable a world of 
fascinating experiments in nonlinear 
dynamics. ) 
e A supplemental reference guide 
which is substantially more forthcom- 
ing about the details of the numerical 
algorithms. As I have insisted, informa- 


352 COMPUTERS IN PHYSICS, MAY/JUN 1991 


providing pulley capabilities—a more 
difficult proposition. ) 

e Numerically specifiable surface and 
elasticity parameters. (These omis- 
sions—the result of a conscious deci- 
sion to shield users from the quantita- 
tive aspects of models deemed too 
unphysical to merit the costume of 
numerical precision—only diminish 
one’s ability to use the program intelli- 
gently and reproducibly.) 

@ A display of rotational inertia in the 
spec sheet of mass objects, and the 
option of adjusting it up to the physical 
limits imposed by the mass and physical 
size of the object. (At present, the 
rotational inertia about its center-of- 
mass is calculated correctly—even for 
arbitrarily shaped polygons—under the 
assumption that the object’s mass is 
uniformly distributed over its area; 
however, it is not displayed or separate- 
ly adjustable. ) 

e@ Velocities specifiable by magnitude 
and direction. (For instance, in projec- 
tile motion experiments, one would 
really like to change the initial angle of 
projection without altering the speed.) 
e Meters which are able to display 
vector quantities in magnitude/direc- 
tion format. 

e@ Keyboard equivalents for virtually all 
menu items, as well as for selecting the 
drawing tools. 

e A “Skip Frames” option that is active 
during the initial run. (At present, 
animation frames can be skipped only 
during replay. The speed and accuracy 
of some simulations would benefit dra- 
matically if one could reduce the inte- 


gration time step and request that IP 
refrain from drawing most frames.) 
e All sorts of meters beyond the 
admittedly substantial number already 
provided. In particular, I would like a 
“Total Energy” meter which could be 
used to monitor experiments devoid of 
friction and inelasticity as a crude check 
for the build-up of numerical errors. 
è Optional numerical control over the 
positions of the no-mass objects. The 
user should be able to specify explicit 
positions for the “business ends” of all 
no-mass objects. (In the case of strings, 
springs, and dashpots, one can at least 
force them to be connected to objects at 
their centers-of-mass; no such facility is 
provided for force objects.) 

è Reshapeable polygons. 

Next, a wish list of items which 
would most likely take longer to be 
implemented: 

è Pulleys. (This is not as simple as it 
may sound since it would require 
strings to be imbued with the ability to 
detect contact with mass objects along 
their entire periphery. As soon as we 
begin thinking along these lines, we will 
want nonmassless strings which sup- 
port waves. Lots of calculations going 
on now!) 

è Sticky surfaces that can support an 
attractive normal force, perhaps up to 
some maximum value. (This is also 
likely to be harder than it sounds, due 
to the complications it will introduce 
into the friction and elasticity algor- 
ithms. ) 

e A much better air resistance algo- 
rithm. Air resistance is an interaction 
that is exceedingly difficult to model 
correctly. To do so and maintain any 
real-time animation capability is nearly 
impossible. JP makes some serious 
compromises here and implements air 
resistance using an algorithm in which 
the force is given by F,,, = — kvl? 
where k is the air resistance parameter 
set by the user (and applicable to all 
objects in any given simulation), v is 
the translational velocity of the body’s 
center-of-mass, and 11 is the length of 
the body’s projected section perpendic- 
ular to its direction of motion. The force 
is directed opposite the velocity and is 
applied at the center-of-mass. 

As a result, semi-realistic results 
are possible only in the case of non- 


rotating bodies which possess and 
maintain symmetrical aspects in the 
direction of motion. For instance, a 
rotating object will experience no decel- 
erating torque and a moving asymmet- 
rical object will neither be deflected nor 
rotated by the effects of air resistance. 
Don’t even think of trying to simulate a 
knuckle-ball! 

On the other hand, the algorithm 
is adequate enough for simulating the 
motion of a simple projectile or the 
deceleration of a disk moving along a 
frictionless surface. Even accepting its 
current limitations, the manual should 


e A gravitational interaction between 
objects, perhaps with variable depen- 
dence on distance, for exploring orbital 
mechanics. (Here I must briefly men- 
tion Gravitation, LTD., a remarkable 
shareware program written by Jeff 
Romereide, 343 Elma Avenue, Laurel 
Springs, NJ 08021. Although it deals 
with a more limited realm of phenome- 
na than IP, Gravitation, LTD. is an- 
other must-have piece of software.) 
The Knowledge Revolution people 
are well aware of all these desires and, 
not surprisingly, have dreams of their 
own that go well beyond. Nevertheless, 


IP can also become the research 


. instrument for an unlimited number of 


student-guided mini-research projects. 
For instance, Fig. 1 shows a frame from 
a quickly designed and run experiment 
to investigate the statistical behaviors 
leading to the equipartition of energy in 
an ideal gas. A “gas” consisting of eight 
1-kg disks and one 4-kg disk is confined 
inside a container. The initial velocities 
of all objects except the 4-kg disk are 
zero. The 4-kg disk is given an initial 
velocity of 5 m/s (KE = 50J) toward 
the others. All collisions are frictionless 
and perfectly elastic. Running the ex- 





Fig. 1. A frame from an experiment designed to explore equipartition of energy. Eight 1-kg disks are initially 
at rest, but confined within a gravity-free box with a single 4-kg disk which has an initial speed of 5 m/s. All 
interactions are perfectly elastic and frictionless. (The total kinetic energy of the system is 50 J or about 5.6 J 
per object.) 


disclose the algorithm being used, its 
severe limitations, and the units of k 
(which are kg/m’). 

@ A way in which to make the gravita- 
tional field vector point into the screen 
with objects sliding and rotating fric- 
tionally on a two-dimensional surface. 
This is a highly non-trivial problem. 
The frictional forces and torques exert- 
ed by a surface on a rotating and 
translating body are very complicated 
as a result of the fact that each element 
of the contact surface is moving in a 
different direction. Of course, this is 
why we want a computer simulation to 
do it for us! 

@ The ability to give objects arbitrary 
charges—or, better yet, charge distribu- 
tions—and to specify external electric 
and magnetic fields. As long as we’re at 
it, let’s have a magnetic dipole moment 
as well. 


it is important to keep in mind that 
substantial new capabilities will almost 
certainly come at the cost of both 
simple operation and animation speed. 


What Is It Good For? 


First and foremost, JP deserves a place 
in every mechanic’s classroom as a 
device that bridges the gap between the 
blackboard and traditional demonstra- 
tions. The blackboard is a static medi- 
um, woefully inadequate when it comes 
to conveying dynamic information. On 
the other hand, demonstrations are 
notoriously treacherous, and even on a 
good day it’s not possible to vary all the 
parameters one might like. With a 
projection system and a moderate in- 
vestment of time, an instructor can use 
IP to great advantage as an in-class 
supplement. 


periment provides a dramatic, visually- 
delivered lesson on how a gas is “‘ther- 
malized.” It argues poignantly in favor 
of the second law of thermodynamics, 
and it makes clear that the average 
speed of the large disk and the relative 
variations of its speed are smaller than 
those of the others. 

Furthermore, by metering the 
large disk and one of the small disks, 
and exporting the data for external 
analysis, one can perform calculations 
which augment the observations. Figs. 
2 and 3 show, respectively, the speeds 
and kinetic energies of the two disks as 
functions of time and the averages and 
standard deviations of these quantities. 
Note that the average speed of the small 
disk is twice that of the large one, while 
the energies are more comparable. The 
average energy of the large disk, which 
experiences smaller variations, is very 
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PRODUCT REVIEW 





Particle Speeds 
— 1kg(3.3+ 2.0) m/s 
women 4 kg (1.6+0.6) m/s 





Time 


Fig. 2. Data from the experiment of Fig. 1. The speeds of one of the 1-kg objects and the 4-kg object 
are plotted as functions of time. The legend gives the average and the standard deviation of the val- 
ues for each object. The first few data points have been left out of the statistical calculations to allow 
the system to “thermalize.” This figure and Fig. 3 were created by exporting data from /P into 
Cricket Graph.” 


30 


Particle Energies 


— 1kg(7.4+7.6)J 





Time 


Fig. 3. Data from the experiment of Fig. 1. The energies of one of the 1-kg objects and the 4-kg 
object are plotted as functions of time. (See the caption for Fig. 2.) 


close to one-ninth the total kinetic absolute and relative variations in the 
energy of the gas (50 J/9 = 5.6 J). A speeds and energies. 

number of other statistical predictions Computer-based research projects 
can be checked by looking at the like this combine the essential elements 
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of experimental and theoretical re- 
search. They provide time-efficient 
methods for students to explore phe- 
nomena and to collect data unencum- 
bered (for the moment) by enormous 
instrumental roadblocks. The data can 
be analyzed empirically for trends 
and/or examined for conformity with 
theory. Although the example of the 
last paragraph was intended to show off 
some of the abilities of JP, and might be 
better suited to students who have 
completed their introductory physics 
courses, similar but less ambitious pro- 
jects are appropriate even for beginning 
students and can foster an early appre- 
ciation of the excitement of indepen- 
dent research. 


What Do You Think Overall? 


All things considered, I like ZP. Al- 
though my wish list is long, my admira- 
tion for the capabilities already there is 
substantial. Yes, it has bugs and, yes, 
even beyond the bugs it can easily get in 
over its head. One has to be aware of its 
limitations and know how to work 
around them. On the other hand, more 
than once have I learned some physics 
by spending a little time thinking about 
an odd result that I had believed at first 
to be the result of a flaw in the program. 
And even when it is wrong, valuable 
lessons can be learned from considering 
the cause of the failure; there are 
remarkable subtleties involved in apply- 
ing Newton’s laws to the interactions of 
extended bodies. The more I use JP, the 
more things I think of to do with it. It 
supports a wide range of applications, 
and I can imagine no readier form of 
simulation set-up. Get your hands on it 
and come to your own conclusions, but 
if you are an instructor of high school or 
lower-division college physics—or if 
you simply enjoy a headier sort of 
computer game—you will want to learn 
about this program. a 
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